Cusp Points in the Parameter Space of 
Degenerate 3-RPR Planar Parallel Manipulators 
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This paper investigates the conditions in the design param- 
eter space for the existence and distribution of the cusp lo- 
cus for planar parallel manipulators. Cusp points make pos- 
sible non-singular assembly-mode changing motion, which 
increases the maximum singularity-free workspace. An ac- 
curate algorithm for the determination is proposed amend- 
ing some imprecisions done by previous existing algorithms. 
This is combined with methods of Cylindric Algebraic De- 
composition, Grobner bases and Discriminant Varieties in 
order to partition the parameter space into cells with con- 
stant number of cusp points. These algorithms will allow us 
to classify a family of degenerate 3-RPR manipulators. 

Keywords: kinematics, parallel manipulator, singulari- 
ties, cusp, discriminant variety, cylindric algebraic decom- 
position, degenerate 3-RPR, symbolic computation. 



1 Introduction 

In the past, singularities were believed to physically sep- 
arate the different assembly modes, meaning that for fixed 
joint values one could not find a path going from one assem- 
bly mode to another without crossing a singular configura- 
tion. So the interest relied on considering the widest con- 



nected non-singular domain, called aspect. Innocenti and 
Parenti-Castelli pointed out in [ 1 1 that non-singular changes 
of assembly mode are possible, and McAree and Daniel 
showed in fT] that such changes are possible when triple 
roots of the Forward Kinematic Problem (FKP) exist. In O 
Zein, Wenger and Chablat showed that for the case of 3-RPR 
manipulators a non-singular change of assembly mode can 
be accomplished by encircling a cusp point, and Husty re- 
cently proved in J?) that the generic 3-RPR parallel manipu- 
lators without joint limits always have 2 aspects. 

From the algebraic point of view, the locus of cusp 
points can be described by means of symbolic equations. In 
order to avoid long symbolic-algebraic manipulations, these 
equations are usually solved by numerical approximation at 
an early stage, which may lead to small deviations that can 
be propagated along the process. However, there exist effi- 
cient symbolic-algebraic techniques that may leave the use 
of numerical methods to the last step. In particular, we will 
apply Grobner bases |5| in order to adopt a more suitable 
equivalent system defining the same solution points. 

Lazard and Rouillier introduced the mathematical no- 
tion of Discriminant Variety (DV) |j6l, which is a variety of 
codimension 1 in the chosen parameter space whose com- 



plement satisfies the property that over each connected com- 
ponent the given system has a constant number of solutions. 
The complement of this DV will be partitioned into cells by a 
Cylindric Algebraic Decomposition |7|, also known as CAD. 

This paper is intended to illustrate both the performance 
of the new algorithm for the determination of the locus of 
cusp points and its combination with the forementioned al- 
gebraic techniques in the analysis of existence conditions 
and distribution along a 2-dimensional parameter space. Al- 
though the method can be applied to more general manip- 
ulators (see [8 1), such performance will be exemplified on 
a family of degenerate 3-RPR manipulators, detailed in Sec- 
tion|2] The algorithm for the cusp point determination, which 
is one of the main contributions of the paper, is given in Sec- 
tion [51 where it is compared to other previous algorithms. 
Section |4] outlines some of the exploited algebraic objects 
such as the DV. In Section |5] the previous procedures are 
combined with the CAD to partition a 2-dimensional space 
with regard to the associated number of cusp points, which 
leads us to analyze a complete family of degenerate 3-RPR 
manipulators that depend on one geometric parameter This 
section also illustrates some applications of the presented 
strategy to robot design. The paper concludes in Section |6] 

2 A class of degenerate 3-RPR 

Let us describe the family of manipulators on which the 
strategies presented along the paper will be exemplified. A 
general 3-RPR manipulator is a 3-degrees-of-freedom pla- 
nar parallel mechanism that has two platforms connected by 
three RPR rods, with the prismatic joints being actuated and 
the revolute ones being passive. Without loss of generality 
we can assume the absolute reference frame to be such that 
the base points of the leg rods are Ai = (0,0), A2 = {A2.x,0) 
with Ajx > 0, and A3 = (A3^,A3y). If Bi, B2 and B3 are the 
corresponding points on the moving platform, then the geo- 
metric parameters associated to this manipulator are the val- 
ues A2x,A3.:c,A3v, the lengths di = ||BiB2||, "^3 = ||BiB3||, and 
the angle p = B2B\Bt,. The input-space is then formed by 
p = (pi,p2,P3) € M^, where p,- > are the leg rod lengths, 
and the output-space is formed by the poses of the moving 
platform x = (x,}', a), where B\ = {x,y) and a is the angle 
of vector B2 — By relative to A2 —A\. We define (ia,Ca), 
(ipjCp) and (^a+pjCa+p) to denote the sines and cosines of 
a, P, and (a + P), respectively. Then the forward kinematics 
of a general 3-RPR manipulator is defined by the system of 
equations 

{x + diCa-A2^f + {y + diSaf-9l=Q (1) 
{x + di Ca+p - A3x)^ + {y + di .Sa+p - A3jO^ - P3 = 0. 

For these manipulators. Hunt showed that the FKP ad- 
mits at most 6 assembly modes f9l, and several authors f TO] 
QlJ proved independently that the system associated to the 
FKP can be reduced to a polynomial of degree 6. The 3- 
RPR manipulators for which the degree of this characteris- 




Fig. 1 . Example of degenerate 3-RPR. 



tic polynomial decreases are known as analytic or degener- 
ate I 12I13II . because the Cramer system in Gosselin's method 
degenerates. In this paper we will focus on a class of degen- 
erate 3-RPR manipulators whose base and moving platforms 
are congruent triangles, with the moving triangle being re- 
flected with respect to the base one, as that of Fig. [T] This 
class of manipulators was first studied by Wenger, Chablat 
and Zein in |14|. Their mathematical description requires 
the addition to the initial Eqn. ([T]) of the following geometric 
constraints: 

dl = A2;t 

cos(p) = Ai,/di (2) 
sin(p) = -Ai,/d3. 

Therefore, the system of equations defining this family of de- 
generate 3-RPR manipulators, formed by Eqns. ^ and (|2), 
will be denoted as 

F(p,x)=0. 

Generically, we will refer to this system by F . Whenever 
concrete values for the geometric parameters are considered, 
these will be specified. Finally, the notation will stand 
for the evaluation on real values (p,x). 

3 Cusp locus determination 

In this section we describe the cuspidal locus and an- 
alyze the usual methods for their determination. After that 
we propose a more accurate approach and compare it to the 
previous ones in a simple example. 

Let us assume that a specific manipulator, whose geo- 
metric parameters have been set into F, has been designated. 
Then, we denote the associated configuration space by 

ai^) = {(p,x)eM^:F|(p,,)=0}. 

The Jacobian matrix of F with respect to the output vari- 

/ -jr- "lE- liE' \ 

ables is denoted as Jx(^') = ( "gj "5" gj;' ) ■ The configura- 



tions where its determinant is zero are called parallel singu- 
lar configurations, or type 2 singularities. On these configu- 
rations the manipulator shows a loss of control. The parallel 
singular locus of our manipulator is a 2-dimensional space 
that can be described (see 1 15 1) as 

1(F) = {(p,x) e C{,F) : Jx(^')|(p,x) is rank deficient}. 

For simplicity, we will refer to this set as the singular locus. 
With this setting we now define the cuspidal locus as 

k(F) = {(p,x) e C{F) : p root of exact multipHcity 3 of F}, 

i.e. the tiiple roots of the FKP. Observe that k(F) C E(f ), 
since the Jacobian Jx(^^) is rank deficient on the roots of mul- 
tiplicity three of F . It is known that in the proximity of cusp 
points a non-singular change of assembly mode can be made. 
Figure |2] shows a cusp point K and a non-singular path con- 
necting two different assembly modes (pi and p^). We shall 
note, however, that both the singular and the cusp locus are 
quite difficult to visualize in the 6-dimensional (p,x)-space. 
So for mechanisms with at most one inverse kinematics so- 
lution, as is the case for the 3-RPR, we will actually project 
them onto the input-space instead. 

3.1 Usual methods for the cusp computation 

Let us revise the two main algorithms that have been 
more commonly used in the determination of the locus of 
cusp points for a given manipulator The following method, 
introduced by Wenger and Chablat in ifOl . and anallytically 
derived in 1 16 |, has been used for the degenerate 3-RPR ma- 
nipulators. It was inspired on an approach developed by 
Hernandez et al in ifTTl for other robots. 



Algorithm 1 by Wenger and Chablat |T3l 

1 . Reduce F (by successive resultants) to a single equation 

g{t) = 0, with t = tan(a/2) and coefficients in p. 

2. Equations of triple roots of g 

G-U = 0,|=0,||=0}. 

3. Equations of strictly triple roots of g 

G = GU{|f ^0}. 

4. G — Eliminate t from G and solve the remaining system 
for real values of p. 

5. Solve G. 



This strategy reduces the problem to the computation of 
the strictly triple roots of one single univariate polynomial g. 
However, the constraint added in step 3 makes the computa- 
tion quite hard, and thus this step is often removed. 

Another commonly used method, described by McAree 
and Daniel in 121, makes use of the series expansion of F. 




Fig. 2. Cusp point K as a triple root of the FKP and non-singular 
path linking upper and lower solutions of the FKP. 



Algorithm 2 by McAree and Daniel |2| 

1 . Series expansion of F 

= f Ax + |Ap + • Ax^ (g) Ax + Ax^ Ap+ 

|V(0)ap + ... 

2. Compute configurations where 1st and 2nd order con- 
straints are rank deficient, i.e. solve 

^u|^^ V = 0, where v is a unit vector in right kernel 

of and u is a unit vector that spans left kernel. 



This second strategy reduces the problem to the resolu- 
tion of some quadratic equations, but it also requires to find 
the unit vectors u and v, which may hinder the computation. 

These algorithms are commonly used in the cusp locus 
determination. However, both have drawbacks related to the 
non-cuspidality of some resulting points: 

• Since the polynomial g obtained by Algorithm [T] is the 
result of several projections, some of the obtained points 
may correspond to the projection of complex (not real) 
solutions, as we will see later on. 

• Step 3 usually needs to be removed from Algorithm [Tjin 
order to avoid slow-processing. 

• Algorithm |2] does not constrain the multiplicity of the 
solutions to be exactly 3, so it may obtain higher multi- 
plicity ones. Regardless of that, in 1 3 1 it is shown that ad- 
ditional spurious solutions may be produced for generic 
3-RPR manipulators. 

Therefore, both methods can only provide sufficient condi- 
tions for the cuspidal locus but not always necessary ones. 

3.2 Improved method 

Despite the fact that the formulation of the cusp locus 
is quite simple, the associated system of equations usually 
contains many equations in many unknowns, whose resolu- 
tion can take long computations and even lead to abnormal 
termination for not too complex examples. So the methods 
described previously were introduced as simple, though not 
accurate, alternatives to the symbolic resolution. However, 
we can now get over some of these difficulties with current 
powerful symbolic algebra tools that fix the deficiencies of 
the algorithms detailed above. 



The approach that we propose is an evolution of ifTSl 
by Moroz et al., inspired on the results of |19|. The main 
difference of the proposed method compared to that of ifTSi 
is the introduction of the saturation operator to remove the 
quadruple roots. 



simple case of degenerate 3-RPR in order to contrast their 
results. However, let us clarify that both the formulation 
and the proposed algorithm apply to other more general ma- 
nipulators (see lis]). We set the geometric parameter values 
Aix = l,A3x = 0, A3v = 1, P = -7t/2, di = 1, andc/3 = 1. 
The characteristic polynomial for Algorithm [Tjis 



Algorithm 3 Proposed method 

1. Equations of double roots of F w.rt. p 
DF=FU{det{UF))=0} 

2. Equations of triple roots of F w.rt. p 

7> = Dp U {det(m) : m maximal minors of Jx (£>/?)} 

3. Equations of quadruple roots of F w.rt. p 

Qp — TpU {det(m) : m maximal minors of Jx(7/?)} 

4. Saturate 7> by gf 
CF = satiTF,QF) 

5. Solve Cf for real values of (p,x) 



g{t) = (p2-p2)f3 + (p2_p2_4)^2^(p2_p2_4)^^p2_p2 

After eUminating t from G we get G = {Pi ,^2,^*3} as follows 

p^ + pt-2p2p2+6p2-3p^-3p^-12 = 
2pt + 2p4-4p2p2+4p2 + 3p^-7p2-16 = (4) 
p3 + P?P2-P?P3-p2p3+3p? + pi-4pi-6 = 0. 

Observe that these equations are not independent. Indeed, 
is a combination of the other two: 



Given the system defining the mechanism F, it computes 
iteratively the equations Tf and Qf of triple and quadruple 
roots p of F, respectively. Then, we use saturation. Given 
two polynomial systems and ^2, sat(5'i , 52) is an algebraic 
operator that returns a polynomial system whose solution set 
is the closure of the solutions of the first system after remov- 
ing those of the second one. If y(5', ) denotes the solution set 
of Si, it is satisfied that 



y(sat(5i,52))=y(5i)\y(52). 



(3) 



In general, the saturation ensures that all roots of S2 are re- 
moved. However, in specific cases, some points can remain 
due to property of Eqn. ^ for which we can only obtain 
y(5i)\y (52) instead of y (5i ) \ y (^2), which can differ by a 
null-measure set that can easily be removed afterwards. Fur- 
ther details on the saturation and its geometric interpretation 
can be found in |l5l. 

Although we are only interested in real (feasible) solu- 
tions, we shall note that the polynomial system obtained after 
saturating has real coefficients and thus its solution set could 
contain some complex (not real) roots. For this reason we 
need to solve the final cusp system Cf in the real field. This 
is done by using the RootFinding Maple package. 

With Algorithm|3]the previous drawbacks are amended: 

• When computing the saturation of Tf by Qf, the points 
of multiplicity 4 or higher are removed, and so we can 
guarantee that only the cusp locus is obtained. 

• By solving Cf for (p,x), instead projecting onto the p- 
space, we avoid having biased points produced by the 
projection of complex (not real) solutions. 

• Furthermore, solving Cf in the real field ensures that no 
other spurious complex solutions are considered. 

3.3 Case study comparison 

Let us now compare the performance of both Algo- 
rithm [U without step 3 and the proposed Algorithm |3] on a 
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Additionally, there are solutions of G that do not correspond 
to the real cusp locus. For instance, if we set pi = 1/3, the 
system G|pj = i/3 has two solutions with both p2 and p3 pos- 
itive. But the FKP evaluated on these two solutions only 
has complex solutions {x,y,a). So, for pi = 1/3 the c-space 
C{F) has no cusp points, though Algorithm 1 obtained two 
mistaken candidates. 

We now test Algorithm [3] on the same example. The 
equations of the cusp locus Cf are: 



6ca + P2-P3 =0 

24+%- 1 =0 

2c^-%- 1 =0 
2coria ^ Ca = 
3ca + 3.?a + P?-P3 + l =0 
3ca + 3%+x2+/-p2 + l = 

2sapl-3ca-6sa-4xy-x-y = 
2xca + 2ysa + ca-3sa-2x+l =0 



2yca + 2ysa + ca-sa-x + y+l =0 
4y%-^a + 2(ca-l)p3+4/-33; + x+l=0 
2sa{2y^-4y-l)-3ca-4xy~2y^ + 3y 

-X + P3-2 = 
6sa (2y+l) + 12ca-8/x + 8/ + 12jt;3; 

+x-33' + 2xp3-63'P3-6p3 + 8 =0 
18%(18y- l)+36ca + 32/ + 20/ + 44x3; 
+ 13jt:-473/ + 4p^-8p^yx-32/p^ 

-6xpl -6ypl- 32 p^ + 40 = 0. 



(5) 




Fig. 3. Singular curve for pi = | on (p2,P3) 

Here if we try to solve CF\p^=ip, we get no real solutions. 
Figure |3] shows a section in the (p2,P3)-plane of the singular 
locus of F for pi = ^. As can be observed, there is no cusp 
point in this plot. This phenomenon is not casual. Actually, 
the value 1/3 for p i has not been randomly picked as we will 
see in next section. In fact. Algorithm [3] describes the cusp 
locus more accurately than Algorithm[T] in general. 



4 Discussion on the joint space 

We now extend our improved method to partition a pa- 
rameter space with regard to the associated cusp locus. So, 
we want to discuss the solutions of a parametric system. 
Among the numerous possible ways of solving paramet- 
ric systems, we focus on the use of Discriminant Varieties 
(DV) |6| for two main reasons: it provides a formal decom- 
position of the parameter space through an exactly known 
algebraic variety (no approximation), and it has been suc- 
cessfully used in similar problems [20J. 

Let us consider a general parametric polynomial system 

5 = {pi (v) = 0, . . . ,/7„,(v) - 0,<?i (v) > 0, . . . ,^z(v) > 0}, 

where pi, . . . ,p„,,qi, . . . ,qi aie polynomials with rational co- 
efficients depending on v — {Ui, . . . ,Ud,Xi, . . . ,X„) with X, 
being unknowns and [/, parameters. For instance, the system 
describing the cuspidal configurations our manipulator Cp is 
parametric if some of the geometric parameters are initially 
left free in F. The DV associated to system is described 
by a polynomial equation. This DV partitions the parameter 
space into several regions such that over each open region 
delimited by the DV the number of real solutions of is 
constant. Prior to defining the DV associated to !f, we need 
to specify a solver of 0-dimensional systems that will be used 
as a black box. 

4.1 Basic black-boxes 

Let us describe the global solver for 0-dimensional sys- 
tems that will be used as a black box in the general algorithm. 
We mainly use exact computations, namely formal elimina- 
tion of variables (resultants, Grobner bases) and resolution 
of 0-dimensional systems, including univariate polynomials. 



We first compute a Grobner basis of the ideal 
{pi, . . . , p,„) for any ordering, which will help us detect if 
the system has or has not finitely many complex solutions. If 
yes, then compute a so called Rational Univariate Represen- 
tation (RUR) of {pi,... ,pm) (see 121]), which is an equiva- 
lent system of the form 

{.f{T) = 0,^1 = 7(Tf ' • • • ~ g^TT'J'' 

where T is a new variable independent of Xi,...,X„, 
equipped with a so called separating element (injective on 
the solutions of the system) m e Q[Xi , . . . ,X„] and such that : 

y(/7i,...,/7,„) A y(/) ^ v{pi,...,p,n) 

{xi,...,x„) ^^ = u{xi,...,x„) ^ (|(P),...,|(P)) 

defines a bijection between the (real) roots of the system and 
the (real) roots of the univariate polynomial /. 

We then solve / = 0, computing so called isolating in- 
tervals for its real roots, i.e. non-overlapping intervals with 
rational bounds that contain a unique real root of / (see [22l)- 
Finally, interval arithmetic is used in order to get isolating 
boxes of the real roots of the system (non-overlapping prod- 
ucts of intervals with rational bounds containing a unique 
real root of the system), by studying the RUR over the iso- 
lating intervals of /. 

In practice, we use the function RootFinding[Isolate] 
from Maple software, which performs exactly the compu- 
tations described above. 

4.2 Discriminant varieties 

Consider now the constructible set 

5 = {veC":/9i(v)=0,...,/9,„(v)=0, 
^i(v)^0,...,^,(v)^0}, 

and let us assume that for almost all the parameter values 
this 5 is a finite set of points. Then, a discriminant variety of 
S with respect to {Ui,...,Ud) is a variety 1^ C C' such that 
over each connected open set Zl not intersecting 1/ (UO'U — 
0), S defines an analytic covering. In particular, the number 
of points of S over any point of Zl is constant. 

Discriminant varieties can be computed using basic and 
well-known tools from computer algebra such as Grobner 
bases 16J. A full package is available in Maple software 
through the RootFinding[Parametric] package, which pro- 
vides us with a polynomial DV{S', Ui,..., Uj) whose associ- 
ated discriminant variety is 

4.3 Case study comparison 

Let us consider again the degenerate 3-RPR with the 
same geometric parameter values as those specified in sec- 
tion (33] i.e. = 1, A3., = 0, A3,. = 1, p = -7t/2, £/i = 1, 
and ii3 = 1, and consider the systems G (Eqn. |4|, and Cf 
(Eqn.|5]), obtained by Algorithm [T] without step 3 and by Al- 
gorithmic] respectively. We will regard as a parameter one 
of the leg lengths pi of the manipulator The discriminant 
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variety will provide us with a polynomial in p i whose roots 
will delimit some open intervals such that for whatever value 
of pi within one interval, the number of cusp points of F|p, 
wiU be the same. We compute the DV for each system with 
respect to pi, and analyze the results. 

In the first case, we get the polynomial DV{G;^i) — 
Pi (pj — 2) (2p[ + 10p[ — 1), whose roots describing the dis- 
criminant variety are 



ro=0. 



n 



^2^27- 10 
2 



and = ^/2. 



Since the number of cusp points is kept constant between two 
consecutive roots, we can compute the associated number of 
cusps by picking one single value of pi inside each open 
interval and solve Gjp, . In this case we obtain 

• cusp configurations for pi e ]0, ri [, 

• 2 cusp configurations for p i € ] r i , rs [, and 

• 3 for pi e]r3,oo[. 

Substituting p i = r, into G we obtain the number of cusps on 
the borders of the intervals. 

• cusps on pi =0, 

• 2 cusps on p 1 = r\, and 

• 3 on pi = r3. 

In the second case, a similar analysis for Cf gives 



DV{Cf\9i 



Pi(Pi 



2)(8p?- 



l)(2pt- 



10p?-l), 



which has one more root than DV (G; p i 



ro = 0, ri = ^ — 2 ' 



r2 = ^ and 7-3 ^ v^. 



• cusps on pi =0, 

• cusps on pi = r\, 

• 2 cusps on pi = r2, 

• 5 cusps on pi = r3. 



The intervals and the numbers of cusps for Cf differ a bit 
from those obtained for G: 

• cusps for pi e ]0, r\ [, 

• cusps for pi e]ri,r2[, 

• 4 cusps for pi e]r2,r3[, 

• 6 cusps for pi e ]r3,oo[. 

The results obtained in both cases are compared in 
Fig. |4] We can observe that the first two intervals do not 
exactly coincide, and that for the second system the obtained 
numbers of cusps appear doubled for all intervals (compared 
to those obtained for the first system). Both phenomena can 
be explained as a consequence of the projection map used to 
compute the system G. Let us remind the reader that G is 



obtained after several reductions of the initial system, each 
of which applying also a projection on the p-space. For this, 
there can be complex configurations of the manipulator that 
project onto real roots p of G. This is the case for the values 
of pi e]ri,r2[. 

The same can be done for any other parameter and the 
same phenomena can be observed. 



5 Higher-dimensional discussion by means of a CAD 

By construction we know that over any connected open 
region not intersecting the DV the system has a constant 
number of real roots, for whatever chosen parameters. But 
if we want to discuss larger parameter spaces, then the open 
regions will no longer be as simple as 1 -dimensional inter- 
vals. So the goal of this section is to provide an accurate 
description of the regions with constant number of solutions. 
For this we will use the Cylindric Algebraic Decomposition 
(CAD) I 



5.1 The complementary of a discriminant variety 

Let C Q[t/i , . . . jf/f/] be the set of polynomials de- 
scribing the DV. Then for each / = c/ — 1 , . . . , 0, we introduce 
a new set of polynomials C (Q)[t/i , . . . , Ud-i] defined by a 
backward recursion: 

• = the polynomials defining the DV, 

• 'Pi = {DV{p;Ui), LeadingCoefficient(/7, U,), 

Resultant(/5,g',f/,), p,q^ !P.^,} 

Each Ti has an associated algebraic variety of dimension 
at most i — 1, "Vi — Vdlperp, /')- The 1^ are used to recur- 
sively define a finite union of simply connected open subsets 
U^!^[ Ui^k C M' of dimension / such that V,- n = 0. 

Before defining the sets llik, we introduce some nota- 
tion: for a univariate polynomial p with n real roots. 



root(p,Z) 



-c=oif / < 0, 

the real root of p if 1 < Z < «, 

-00 if / > n . 



Moreover, if p is a n-variate polynomial, and v is a (n — 1)- 
tuple, then p^ denotes the univariate polynomial where the 
first n — 1 variables have been replaced by v. 

The recursive process defining the Uij^ is the following: 

• For/= l,letpi =npG!Pi/'- 

Taking all lii =]root(pi,A:);root(/:ii,A: + 1)[ for k = 
0, . . . ,n, where n is the number of real roots of pi, one 
gets a partition of R that fits the above definition. More- 
over, one can arbitrarily chose one rational point mi ^ in 
each open interval Ui,k- 

• Then, for i — 2, ... ,d, let pi — YlpeVj P- 

The regions Zlj ^ and the points Uj ^ are of the form: 

Ui^k ={(vi,...,v,-_i,v,') I V := (vi,...,v;_i) e 

v,e]root(p^O,root(p J, / + !)[} 
Ui.k = (Pi,...,p,-i,P,), with 

(Pl,...,P,_l) =M,-1,; 

p,e]root(p;:"-'^/),root(p:"-'^/ + l)[, 




Pi 

Fig. 5. Plot of the DV of Cp with respect to (p i , t/i ) 
where j,l are fixed integers. 

With this recursive procedure we get a full description 
of the complementary of the DV for the system to be solved: 
the cells IZ^/ and a test point uj i^ £ llj i^ (with rational co- 
ordinates). The number of solutions associated to each open 
cell lld^k is obtained by solving the given system restricted to 
Mrfjt using a 0-dimensional solver Both the cell decomposi- 
tion and the test points can be obtained by the Maple function 
RootFinding[Parametric ][CellDecomposition ]. 

5.2 Open CAD for a class of degenerate 3-RPR 

Let us see the performance of this CAD on a 2- 
dimensional discussion. We consider now a family of degen- 
erate 3-RPR manipulators with A2x = 1, = 0, As,. = 1, 
P = —k/2, and dj, = 1, and regard as parameters both pi and 
di, constrained by di > 0. Now, the system Cf describing 
the cusp locus associated to this family of manipulators has 
18 polynomials, its DV is plot in Fig.|5] and the polynomial 
DV{CF',Pi,di) factors as follows: 

i/iPi(t/? + l)(-4pt-12pt + 27p?t/^ + 15p?-4) 
(4pf + 12pV?-15pi4+4t/f-27p?t/^) 
(256 pt t/? + 8 1 p? c/f - 288 p^ d^ + 256p\- 576 p^ d^i 
+ 51p?4-16^/f-288pt + 51pit/? + 81p?). 

The complement of this DV produces 90 cells with as- 
sociated numbers of cusps varying among 0, 2, 4 and 6, as 
shown in Fig.|6] where black vertical lines delimit intersec- 
tion points of the DV. Let us notice that although all cells 
must be considered disconnected, it is apparent that cells 
are naturally grouped with regard to their number of cusps. 
Additionally, cells with different number of cusps are exclu- 
sively separated by curves of the discriminant variety. Fur- 
thermore, this distribution is consistent with that obtained for 
(ii = 1 in Section 14.31 as can be seen in Fig. [T] However, 
here the section is divided into many more smaller intervals 
whose borders we cannot apriori ensure to be associated to 
a specific number of cusp points. In fact, let us also observe 
that the cells with 2 cusp points (in blue) degenerate into one 
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Fig. 6. Cell Decomposition for (pi,iii) 
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Pi 

Fig. 7. Zoom in of the cell decomposition for (p 1 , t/i ) . Line t/i = 1 
in white. 



single point for c/i = 1 . So we could claim that the case with 
c/i = 1 is a very special degenerate 3-RPR manipulator. 



5.3 Study of the cusp points on the borders of the CAD 

At this point we can only certify the number of cusp 
points in the open cells. This excludes the cell borders. The 
union of all these borders consists of the DV plus the delim- 
itation of intersection points of the DV. However, by defini- 
tion, changes in the number of solutions can only happen on 
the DV. So, we just need to analyze the DV. 

We could try executing a further iteration of the CAD on 
the DV, but the system to be solved turns out to be too com- 
plex and we cannot obtain any results after long computa- 
tions. It is clear that not all points on the DV will correspond 
to the same number of solutions. But we can expect the num- 
ber of solutions to be preserved along the DV between two 
consecutive auto-intersection points of the DV. And although 
it has not yet been proven, many tests have been run on sev- 
eral examples with random points on the DV and all the re- 
sults confirm what the following conjecture infers. 

Conjecture 1. Given a polynomial system F, let 'U be the 
DV of F w.rt two parameters Ui , U2, and let J? be the set of 
its auto-intersection points. Then, the number of solutions of 
F is constant on each connected component of 

The following algorithm analyzes the to study the num- 
bers of solutions on the DV based on the previous conjecture. 



^1 cusps ^1 1 cusps ^1 2 cusps ^| 3 cusps 

4 cMsps ^1 5 cusps ^1 6 cMsps 




(a) (b) 

Fig. 8. Distribution of cusp points on DV {CF\Pi,di) (a), and 
zoom in view on [0, 1 .5] x [0, 1 .5] (b). 



Algorithm 4 Number of solutions of F on the DV 

■V = vai-iety of DV{F-UuU2) 

SI = {auto-intsersection points of 1/} 

for each connected component IZ, of 1^ \ ^ do 

Pi = random point on Hi 

Compute the number of solutions on IZ, as 

the number of solutions of F\p. 
end for 

for each point ^ G ^ do 

Compute the number of solutions of F\g 
end for 



5.4 Complete analysis and applications 

From the results exposed in the previous subsections and 
by joining all the different pieces together we obtain a com- 
plete partition of the 2-dimensional parameter space. 

The execution of Algorithm |4] on DV{CF',pi,di) pro- 
vides the distribution of cusp points shown in Fig. |8] The 
integral picture of the 2-dimensional distribution is given in 
Fig- in It is interesting to notice that there is a continuity on 
the transitions between cells having the same number of cusp 
points, since their common border inherits that same number 
of cusp points. 

Observe also that this distribution has been obtained 
thanks to the DV associated to the chosen parameters di and 
pi, which depends exclusively on these two parameters. In 
particular, DV (Cf ; p i,di) does not depend on p2 nor p3. This 
tends to be erroneously interpreted as: 

"i/ we pick a (pi, di)-point with associated number of 
cusps k, then whatever the values P2 and P3 may take, 
Cf\ (Pi,p2,p3,i:/i) ^'^■^ ^ solutions'" . 
Instead, it should be read as follows: 

"if we pick a {pi,d\)-point with associated number of cusps 
k and fix these values then, among the reachable configura- 
tions there are k cuspidal ones, i.e. CF\(p^.d^) has k solu- 
tions". However, the number of associated cusp points does 
establish a maximum of cuspidal configurations for any val- 




Pi 

Fig. 9. Complete analysis of the cusp points for (p i , t/i ) 

ues p2 and p3. For example, in yellow regions we can have a 
maximum of 4 cuspidal configurations, but depending on the 
values of p2 or p3 there can even be none. In particular, for 
the red regions there are cusp points for all possible values 
P2 and p3. 

Some applications can be derived which may be inter- 
esting from the designer's point of view: 

• It can be helpful in deciding the most suitable architec- 
ture of the mechanism. Let us assume that we want to 
design a 3-RPR manipulator with some given geometric 
constraints such that for a specific task one of the legs 
has to be blocked to a fixed length pi, but the job re- 
quires a large singularity-free workspace. Therefore, we 
may be interested in finding a range ls.d\ of parameter 
values for which the manipulator is cuspidal. 

• It can also be useful for deciding the most suitable 
ranges of leg lenghts for each possible architecture, 
given a specific task. For instance, let us assume that the 
job is set for a non-cuspidal manipulator with parameter 
values A2x = 1, As;^ = Q,Aj,y = 1, p = -7t/2, and c/3 = 1, 
but it requires the largest possible range of the leg length 
pi. Then, the value di can be optimized with this crite- 
rion. Figure[TO]details both the optimal valule d\ and the 
largest possible range Ap 1 for our problem. 

Let us just notice that in both cases the obtained ranges 
can be of varied topology (open, closed, semi-closed, open 
and closed, connected, or even a union of these types). This 
is due to the combination of both the CAD and the study of 
the cusp locus on the DV. 

6 Conclusions 

This paper has introduced both an efficient method for 
the computation of the cuspidal configurations of a mecha- 
nism, and a reliable algorithm that partitions a given param- 
eter space into open regions with constant number of associ- 
ated cusp points. 

The first one is based on a symbolic-algebraic approach 
able to describe the roots of exact multiplicity 3 and a cer- 
tified numerical algorithm that isolates among them the real 
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Fig. 10. Optimal di and Api for non-cuspidal degenerate 3-RPR. 

(i.e. not complex) ones. This symbolic-numeric approach is 
more efficient than other previously existing methods, which 
mainly relied on the approximation of roots of multiplicity at 
least 3 after reducing the initial system to a simpler one and 
projecting it onto the p-space. 

This new method is combined with some algebraic 
tools such as the discriminant variety (DV) and the cylin- 
dric algebraic decomposition (CAD) in order to analyze a 2- 
dimensional parameter space with respect to the associated 
number of cusp points. This second algorithm provides a 
partition of the parameter space into cells with constant num- 
ber of cusp points, which is certified for whatever values are 
picked inside each open cell but not on their borders. Cell 
borders are further analyzed by Algorithm 5]based on Con- 
jecture 1, which still remains unproved. 

Both algorithms have been applied to the analysis and 
distribution of the cusp locus for a family of degenerate 3- 
RPR manipulators, and some applications to robot design are 
also derived. This does not mean that the two given algo- 
rithms are specially designed for this type of 3-RPR mecha- 
nisms. Indeed, they are suitable for more general examples 
since they do not rely on any ad-hoc formulation. Neverthe- 
less for some examples the obtention of results within rea- 
sonable time may not be feasible yet, since there is an im- 
portant symbolic-algebraic part, and thus the more complex 
the initial system is the harder it will be to compute a parti- 
tion on the parameter space. 
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